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Abstract: We investigate what becomes of the translational zero-mode of a five- 
dimensional domain wall in the presence of gravity, studying the scalar perturbations of a 
thick gravitating domain wall with AdS asymptotics and a well-defined zero-gravity limit. 
Our analysis reveals the presence of a wide resonance which can be seen as a remnant of 
the translational zero-mode present in the domain wall in the absence of gravity and which 
ensures a continuous change of the physical quantities (such as e.g. static potential between 
sources) when the Planck mass is sent to infinity. Provided that the thickness of the wall 
is much smaller than the AdS radius of the space-time, the parameters of this resonance 
do not depend on details of the domain wall's structure, but solely on the geometry of the 
space-time. 



Keywords: 



Field Theories in Higher Dimensions 



Contents 



[IL Introduction |] 

|2L Domain walls in 5D space-times |3| 

|2.1| Non-gravitating case - domain wall in the flat space-time || 

|2.2| Background solution for a self-gravitating domain wall S 



^ Scalar perturbations 

|3.1| General setup 

|3.2| Quantum mechanical analogy || 

|3.3| Nature of the spectrum H 



3.4 Partner system and the determination of the eigenfunctions 10 

^ Static potential between sources on the brane 13 

|5L Periodic source on the brane 17 

Conclus ions |19| 

|A] . /i = solutions of the eigenvalue equation of 7ti 19 

|B|. Connection between the spectra of Tii and Ti^ 21 

Eigenfunctions of the partner operator 7^2 22 

|5| . Determination of the static potential 25 

EL Solutions of the coupled system 26 



1. Introduction 

The proposal that our universe could be a topological defect in a higher-dimensional 
space-time, which is at the basis of the brane-world scenarios, was put forward already 
in the 80's [1,2] as an alternative to Kaluza-Klein type theories [3,4]. Contrarily to the 
Kaluza-Klein type setups, in the brane-world models the extra dimensions can be large and 
even infinite. The key idea is that by trapping the matter fields on the defect (a domain 
wall or a string) one can obtain a low-energy effective theory which looks four-dimensional, 
the extra dimensions becoming visible only in the very high energy experiments. The zero 
modes of the higher-dimensional fields, confined on the defect, can be interpreted as the 
four-dimensional particles. In the setup of Ref. [1] the localization was achieved for two 
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types of fields. Firstly, among the excitations of the scalar field forming the wall there 
was a localized zero-mode, which could be seen as a massless 4D scalar particle. Secondly, 
the model naturally incorporated fermions, allowing trapping of the fermion zero mode via 
Yukawa interaction between the fermion field and the domain wall. 

The interest in higher-dimensional theories of this type increased tremendously few 
years ago, when it was discovered that this kind of picture could in fact assimilate gravity. 

Firstly, in [5,6] a model was proposed where our universe is modeled by a thin 3-brane 
embedded in a higher-dimensional space-time with compact extra dimensions; the influence 
of the brane's tension on the geometry of the bulk is assumed negligible. All the Standard 
Model fields are supposed to be trapped on the brane, only the gravity being able to 
propagate in the bulk. As a consequence, the extra dimensions could be fairly large without 
contradiction with the available small-distance Newton's law measurements (at the time 
when this setup was proposed, even the dimensions of order of 1 mm were not excluded; the 
recent experimental bounds [7] require the extra dimensions to be not larger than 160 /im). 
The spectrum of perturbations in the models of this type has been studied extensively [8- 
11], and it was found to contain Kaluza-Klein modes of gravitons, as well as exotic scalar 
particles, called branons [12]. Branons can be identified as the Goldstone bosons associated 
with the spontaneous breaking of the isometries of the bulk, caused by the presence of the 
brane. Massive branons were advocated as dark matter candidates [13]. 

Soon after the proposal of [5,6], the authors of [14] have built two brane world mod- 
els: firstly a compact one, involving two delta-like 3-branes located at the fixed points 
of an orbifolded AdS§ space-time (RSI), and subsequently a non-compact single brane 
model (RS2), obtained by letting the compactification radius in RSI to be infinite. Con- 
sidering the graviton spectrum in these scenarios, they have shown that it was possible to 
reconcile large (or even infinite) extra dimensions with Newtonian gravity. The effective 
four-dimensional spectrum of both the above models was studied in detail by several au- 
thors (see e.g. [15,16]). It was found that in the RSI model, the effective theory contains a 
massless spin- two tensor field (four-dimensional graviton) , an infinite tower of the massive 
tensor fields (Kaluza-Klein gravitons), and a massless scalar field called radion. When 
sending the size of the extra dimension to infinity, thus passing to the non-compact RS2 
model, the radion can be gauged away and disappears from the effective theory [11,15-17]. 
The spectrum of the RS2 model consists of a localized zero mode and a continuum of 
massive modes. Even if there is no gap separating the zero mode from the continuum, 
the law of gravity on the wall is essentially Newtonian, the corrections yielded by the light 
continuum modes modifying this behavior only at short distances. 

The infinitely thin branes used in [14] may be thought as mere approximations of 
smooth structures, thick branes, obtainable as solutions to coupled gravitational and matter 
equations. Regularized versions of the RS2 model were constructed [18-21] along the same 
lines as in Ref. [1]: smooth gravitating domain walls, generating warped geometry, can 
be produced by a scalar field <I> with a potential spontaneously breaking the <I> — ► — <!> 
symmetry. The spectrum of fluctuations in the regular warped backgrounds has also been 
studied, e.g. in [21] (generalizing Bardeen's formalism of gauge invariant variables) and 
in [22] using the light cone gauge formalism. In this case as well, the only localized particle 
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was the massless spin-two tensor field. 

What may seem surprising in the results concerning the RS2 model and its smooth 
generalizations, is that they do not show a slightest trace of the massless scalar present 
in the setup of Ref . [1] . It seems that in the process of including gravity it has somehow 
"evaporated" . However, it is rather strange that a particle would disappear altogether as 
soon as the gravity is switched on, no matter how weakly. Locally, a gravitating domain wall 
should resemble a non-gravitating one, that is the model of Ref. [1]. Presence of a massless 
particle is bound to have physical consequences - for example we expect it to induce a long- 
range potential between sources located on the brane. The physical quantities (such as, for 
example, static potential between sources) must change in a continuous manner when the 
Planck mass is sent to infinity. There should therefore be some imprint of the translational 
zero mode of [1] on the spectrum of scalar perturbations in the presence of gravity, the 
simplest possibility being that it acquires a finite width or, in other words, becomes a 
resonance - either a standard, massive one [23], or maybe a massless quasilocalized mode, 
of a kind encountered e.g. in the GRS model [24]. 

In this paper we investigate the fate of the scalar zero mode in the presence of gravity, 
considering a thick scalar brane model introduced in [18] which possesses both an AdS§ 
asymptotics and a well defined zero-gravity limit. We use the quadratic action for the 
scalar perturbations derived in [22] to study the interactions of these perturbations with 
sources localized on the brane. Our analysis shows that although there aren't any localized 
modes or even narrow resonances, the light modes of the continuous spectrum ensure 
the continuity of the physical quantities. From the point of view of the four-dimensional 
physics, the behavior of these modes can be interpreted as a wide resonance, whose mass 
and width are entirely determined by the geometry of the space-time and do not depend 
on details of the wall's structure, provided that the wall thickness is much smaller than 
the AdS radius i?AdS of the space-time. The influence of these modes is significant at 
distances much smaller than i?AdSi when the distances R are much larger than i?AdS> the 
potential acquires a, 1/B? form and its amplitude is strongly suppressed. We find similar 
results studying propagation of waves emitted by an oscillating source: for the waves 
with frequencies larger than the inverse AdS radius the standard 1/R dependence of the 
amplitude on the distance to the source is recovered in a large interval of distances. 

The paper is organized as follows: In section ^, we remind the construction of the 
thick domain walls - first in the absence and next in the presence of gravity. Section || 
contains the analysis of the spectrum of the system of scalar perturbations of a gravitating 
domain wall. Sections || and |5| are dedicated to the physical applications: in section |||, we 
determine the static potential between two sources mediated by the scalar perturbations 
and in the section [5] we generalize the study to the case of the moving sources. 

2. Domain walls in 5D space-times 

2.1 Non-gravitating case — domain wall in the flat space-time 

To begin, let us remind in more detail the results of [1], where an explicit field theoretical 
realization of the mechanism of localization of matter on a thick domain wall (without 
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gravity) was devised. The main idea of [1] is to break the translation invariance along the 
extra dimension by forming a domain wall on which the particles will be trapped. This can 
be achieved by a quantum field model containing a real scalar field with the double-well 
potential. The corresponding action is: 



S = d 4 xdr 



(2.1) 



The classical equations of motion admit a domain wall solution $ c (r) interpolating between 
the two classical vacua <£ vac = ±v, whose form coincides with the one-dimensional kink: 

<E> c (r) = utanh(ar) , (2.2) 

where a 2 = Xv 2 /2. The perturbations around this background <f> = $ — <£ c satisfy the 
linearized equation of motion: 

d M d M <j){xi r ) + A(3$ 2 (r) - v 2 )<j){x, r) = (2.3) 
and therefore can be written as 

4>(x,r) = ^u k (x)r] k (r) , 

where u k (x) satisfy the usual four-dimensional Klein-Gordon equation 

[d^ + fi 2 ] u k (x) = 0, 
with mass /i 2 = k v k v and where r\ k (r) are the normal modes of the one-dimensional kink: 

-d 2 r] k (r) + A(3$ 2 (r) _ v 2 )^) = fx 2 Vk {r) . (2.4) 
The spectrum of <fi comprises a zero-mode (fi 2 = 0): 

\/3(2 1 . 

Mx,r) = — u {x) , (2.5) 

I cosh (ar) 

which can be interpreted as a massless scalar particle localized on the wall. There is 
also a localized heavy mode with mass [i 2 = 3a 2 and the continuous spectrum starting 
at fi 2 = 4a 2 , corresponding to the perturbations which are not confined inside the wall. 

For the sake of completeness, let us add that the model can naturally account for the 
massless four-dimensional fermions, localizing the zero mode of the bulk fermion through 
the Yukawa interaction with the domain wall. 

In this work, it is the massless scalar mode (f>o{x, r) and the low energy physics related 



with it which will be the focus of our attention. We expect the particles of the type (2.5) to 
behave as fairly ordinary four dimensional particles: collide, form bound states, etc. Given 
that the particles (j>o(x,r) are massless we expect them to mediate long-range potential 
between sources on the domain wall. 
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2.2 Background solution for a self-gravitating domain wall 

Let us now consider again a theory of a real scalar field $ in a five-dimensional space-time, 
but this time adding an ingredient which was absent in the considerations of Ref. [1]: the 
gravity. Supposing minimal coupling of the scalar to gravity we obtain the action: 



S = id x dr^fg 



R 



1 



~2k + 2 9 



MN 



d M ®d N § - U{$) 



(2.6) 



where k denotes the five-dimensional gravitational constant. In analogy to the zero gravity 
case described in the previous section, we suppose that the potential satisfies U($) = 
{/(—$) and that the field $ is in a domain wall configuration. The presence of the scalar 
field determines the geometry of the space-time via the Einstein equations 



R 



MN 



-gMNR 



i)uW> y i> - um \ ( \g AB d A $d B ® - U($)\ 



and the curvature of the space-time modifies the field equation of the scalar: 

. 



^pd M [y/gg MN d N *]+ Vr 



d<$> 



(2.7) 



(2- 



In order to find the background configuration of the model, one therefore has to solve the 
coupled system of equations Q2.7| ) and Q2.8|). We suppose that the space-time possesses a 
warped geometry, with the metric of the form: 

,2 



ds z 



g M Ndx M dx N = eW^daPda? - dr 2 



(2.9) 



the warp factor e A ^ T ' being a smooth function of the extra coordinate r, with A(r) sym- 
metric and decreasing. The four-dimensional Planck mass is now finite and is defined by 
the relation: 



Mi 



Ml 



dre A ^ , 



where M5 = (4k) 1//3 is the five-dimensional Planck scale. For the metric ( |2.S| ) the back- 
ground equations (|2.7|) and ( pDj| ) read: 

dU 
5$ 



0, 



1/2 



2 k 

y 



.4' 

A' 2 + A" 



12 _ 

2k 

y 



(2.10) 



$' 2 + U 



An interesting exact solution for a thick brane background has been found in Ref. [18]. In 
this solution, the scalar field is supposed to be in the familiar kink configuration: 



3>(r) = v tanh(ar) , 



(2.11) 



which allows to integrate the background equations fl2.10| ), thus obtaining the exponent of 
the warp factor of the form: 



A(r) = —2/3 In cosh 2 (ar) — /3tanh 2 (ar) 



(2.12) 
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and the potential: 



^ A $ 2 ($ 2 



3v 2 



3v 



2\2 



(2.13) 



where again a 2 = \v 2 /2 and where a dimensionless constant (3 = kv 2 /9 has been intro- 
duced. 1 

Besides being an exact solution of ( |2.10| ), the background ( 2.11| )~( |2.13 ) has two at- 
tractive features: 

Firstly, it has a well-defined flat limit, tending smoothly to domain wall solution of [1] 
as (3 goes to zero: 



A(r) 



. 



A 



(<5> 2 -v 



2\2 



It can therefore be considered a self-gravitating extension of [1], which makes it a perfect 
background to investigate the fate of the zero mode of the five-dimensional kink in the 
presence of gravity. 

Secondly, the solution ( 2,11| )-( ^.13| ) can be seen as a regularized version of the RS2 
model [14] and there exists a well-defined "brane" limit (or, in other words, thin wall limit), 
in which it becomes RS2. Indeed, the action ( |2.6| ) can be split into two parts: a localized 
brane part and a bulk part, 



Sbrano = J <^ x dr \/~9 

Stmik = J d * x dr-yfg 



1 



g M »d M $d N $-U($) + U 
2k k 



d 4 xa , 



with the effective bulk cosmological constant A5 determined by the asymptotics of the 
potential (Too = limi,,^^ (7($(r)): 



A5 = kUoo 
and the tension of the brane given by 

dr^fg 

In the thin wall limit, which is defined as [18]: 



2„2 



-24/3^ a 



^{d r ^f + U($)-U 



a — > 00 , 



£ 2 = af3 < 00 



(2.14) 



the warp factor acquires the characteristic AdS form e A ^ r ' — > e _4a/3 ' r and the domain wall 
becomes a delta-like brane, whose tension is tuned to the bulk cosmological constant via 
the relation a 2 = — 6A5/K 2 familiar from the RS models. 



1 Strictly speaking, the potential (2.13) is not bounded from below. However, we are interested in small 
perturbations around the kink configuration and therefore will not be concerned by this instability. 
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3. Scalar perturbations 



Let us now consider perturbations around the background solution ( 2.10 ). When the 
gravity is switched on, we have perturbations coming both from the metric and from the 
scalar field. Their equations of motion can be obtained by setting: 

9M N — ► 9MN + hMN , 

$ — ►$ + </>, 



in the field equations (2/7) and ( |2,8| ) and keeping only terms linear in perturbations. The 
study of the perturbations becomes now fairly involved, as the perturbations of the scalar 
field couple to the metric fluctuations. Yet another difficulty comes from the gauge invari- 
ance - not all the degrees of freedom are actually physical. This inconvenience can be dealt 
with either by constructing gauge invariant variables [21] or by fixing the gauge [21,22]. In 
this work we choose the latter option and use the results of Ref. [22] where the perturbations 
were studied in the light cone gauge. 

3.1 General setup 

In Ref. [22], a formalism was developed to analyze the spectrum of small perturbations of a 
general system of gravitational, gauge and scalar fields in D-dimensional space-times (and, 
in particular, in warped geometries). The bilinear action of the system of perturbations 
was decoupled there into spin-two, -one and -zero sectors using the light cone gauge. In 
the light cone coordinates, defined by x M = (x + , x±, x~, r) = (x + , x 1 , x 2 , x~ , r), with x^ = 
(x° ± x 3 ) / y/2, the background metric ( |2.9[) becomes: 



ds 2 = e A{r) [2dx + dx~ - (dx 



1\2 



(dx 



2\21 



dr' 



The light cone gauge corresponds to setting: 



h. 



M 







for all M 



In this gauge the fields of different spins are readily separated from each other. The physical 
fields are h\2 and h\\ — hyz (gravity multiplet, two degrees of freedom), h ri (vector multiplet 
or graviphoton, two degrees of freedom) and two scalar fields <p and h rr = —h\. 

We will be interested exclusively in the scalar sector of perturbations. It was shown 
in [22] that this sector can be reduced to two coupled scalar fields x{ x i r ) = a/3/8k h rr 
and (p(x, r), where <j) is the fluctuation of the scalar and x a metric perturbation. For D = 5 
their action reads: 



8S. 



(2) 



seal 



/ 



d^xdre 2A ^- 



1 



1/2 



2k& 2 + 



d 2 u \ 

d$> 2 ) 



+^ [e- V'^X " X' 2 + {2A» + A' 2 ) X 2 ] 



.^(2$" + AV)^| • 



(3.1) 
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The fields x an d <j> can be decomposed into four-dimensional plane waves: 

X oc exp(ip^)h p (r) , oc exp(ip M 2^)/ p (r) . 

The four-dimensional scalar spectrum is therefore determined by the following system of 
equations: 



- h" - 2A'ti - (2A" + A' 2 ) h + J— (2$" + & A!) f = ^e~ A h 

V 3 

-/" - 2A'f + + 2A'^ - A") / + ^ (2d>" + VA 1 ) h = ^e~ A f 



(3.2) 



where /j 2 = p^p^ and where we have used the background equation ( p. 10 ) to express 



d 2 u 

in terms of $ and A. For the sake of clarity, we assume that A(r) and $ are given by the 
background solution fl2.11| )-( |2.12| ). 

One should bear in mind that in the light cone formalism the five components of the 
four-dimensional massive graviton are distributed among all three spin sectors: two of them 
are in the tensor sector, another two in the vector sector and the remaining component 
is hidden in the scalar sector. As a consequence, in the scalar sector of perturbations a 
variable can be found which belongs to the graviton multiplet and which therefore has 
the same spectrum as the tensor fluctuations. This variable, which can be written as a 
particular combination of the fields %i 4> an d their derivatives, was identified and described 
in [22]. The main concern of the present work is the physics of the genuine four-dimensional 
scalar field. As the structure of the action makes it clear that it cannot be decoupled using 
a linear transformation, we will work with the coupled system as a whole. 

3.2 Quantum mechanical analogy 

The first observation we can make about the system ( |3.2|) is that in the zero-gravity 
limit {(3 = 0) it decouples and simplifies greatly: 

-h" = ji 2 h , 

which allows us to find its solutions immediately. The solutions of the upper equation are 
simply the plane waves; the second equation, determining /, reduces to Eqn. J2.4Q and 
therefore its solutions are the familiar normal modes of the kink, including the normalized 
zero mode. 

In general, when j3 ^ 0, the problem of finding the solutions of ( |3.2| ) is far less triv- 
ial. Eqn. (|3.2| ) being a system of coupled second order differential equations with fairly 
complicated coefficients, solving it exactly is certainly not an option. It is obvious that 
eliminating one of the functions, passing to a fourth order equation is possible, but by no 
means would it simplify our task. It is preferable to consider the problem of finding the 
solutions of Q3,2| ) as the eigenvalue problem 

W 1 ^ = /x 2 * M , [M , (3.3) 
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for the matrix Schrodinger operator: 



Hi 



dr 2 dr 



\ 



2A" + A' 2 ) 




2$" + §'A') 



- 2A'— + — + 2A' — 
dr 2 dr \ <3?' 



.4' 



\ 



(In the above expressions the index \x denotes the eigenvalue and should not be confused 
for a vector index.) For any given /i, there are four linearly independent solutions of (|3.2|). 
The form of Ji\ implies that its eigenfunctions are normalized with weig ht e A ^ (the precise 
form of the orthogonality relation will be given below). 

The quantum mechanical analogy suggests that we use methods inspired by the su- 
persymmetric quantum mechanics [25]. As we will see, this approach greatly simplifies the 
problem of investigating the spectrum of 7i\. To begin, it can be easily checked that 7i\ 
can be written in a factorized form as 7i\ = A* A, with: 



A = e- A " 



( d 1 ,/ 

dr 4 



V 



At 



-A/4 



( d 1 M 

dr 4 



V 



d L _3 A ,_&^ 
dr 4 ) 

4 <£' / 



= 3A/4 



=3A/4 



where the exponential factors are related with the choice of coordinates and a non-trivial 
weight function present in (EH3). 



3.3 Nature of the spectrum 

Let us investigate whether the spectrum of the operator Tt\ contains any bound states. 
The factorization Hi = A^A guarantees that the eigenvalues of TL\ are non-negative and 
allows one to find easily the /i = solutions of the system (|3.2|) . Let us therefore start 
by checking for presence of a localized zero mode. The zero mode of TCi, if there is one, 
should satisfy either: 

h 
Jo. 



A 



. 



(3.4) 



or an inhomogenous equation 



A 



fo 



(3.5) 



where the inhomogenous term is one of the solutions of A^ ( _° ) =0 (see Appendix Eh . The 

v /o / 

general solution of both the above equations can be calculated exactly. The homogeneous 
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solution reads: 



*o(r) 



h 
Jo. 



( A'e 



(o) 



2k _ , _ i 
3 



+ C, 



(0) 



/ 1-A'e- A 
2k 



dye 



A(v) 



$'e~ 



dye 



A(y) 



(3.6) 



/ 



Let us remark that the solution can be derived as a coordinate transformation: 

5^hMN = —£,N;M — £m-,N 

with = and £ r = e _A . For our background, the solution ( |3,6] ) is clearly not normaliz- 
able; for any choice of cj ^ and we have: 



dre A ^ [^(r)+/ 2 (r)] 



oo . 



This remains true for the two remaining /x = solutions of ( |3.2[ ), which can be determined 
from ( |3.5| ) (as the explicit expressions for these solutions are fairly complicated and not very 
illuminating, we relegate them to Appendix |A], where the determination of the solutions 
of Tii^o = is presented in more detail). There is therefore no localized zero mode in our 
model. 

Let us now turn to the /i / spectrum. To determine the nature of this spectrum, it 
is enough to determine the asymptotic behavior of solutions. Supposing j3 -C 1, far from 
the core of the domain wall, for positive r's, the system (|3.2|) becomes: 



(3.7) 



- h » + 8 fti _ lQ ^ h _ ( fl 2 + <2 ) e -2ar f = ^ 

-/" + 8eY + (4a 2 + 16af) / - 4^ (a 2 + af) e~ 2ar h = /z 2 e^7 , 



where £ 2 = (3a. For |r| 3> 1/a the mixing term is clearly negligible and we are left with 
a simple system of decoupled equations, whose solutions can be written in terms of the 
Bessel functions: 



2e 



o 




+ 



M 2 r 



2? 



where v = The continuous spectrum of 7i\ starts therefore right above zero. There 

are no localized modes. 

3.4 Partner system and the determination of the eigenfunctions 

Absence of bound states in the spectrum of TL\ for (3 ^ does not necessarily mean that the 
localized zero mode present in the non-gravitating case (3 = has altogether disappeared 
as soon as the gravity was switched on. It is very likely that in the presence of gravity the 
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would-be zero mode has become a quasilocalized state. To investigate this possibility, it 
is necessary to know the eigenfunctions of 7i\ for non-zero [i in some more detail. To this 
purpose, it is very convenient to define the "partner" operator TL2 = AA^ . This operator 
is much easier to study, as its eigenvalue problem 




hfj 



has the advantage of being a system of two decoupled equations: 



-h" - 2A'ti 



-A" 



-J" - 2A'f' 



A' 2 I h 



(®'y 5 ... 3 . n 

IF" 2 #0 ~ A V + 2 A + l A 



f 



[j, 2 e~ A h 



» 2 e- A f 



(3.8) 



At the same time, once we know the eigenfunctions of the partner operator H.2, we au- 
tomatically know the eigenfunctions of H\. Indeed, there is a one-to-one correspondence 
between the eigenfunctions of the two operators: for any \i 7^ 0, the solutions of the coupled 
system can be determined from the solutions of the partner system using: 

/ e A l 2 \-h' 




1 

/' 



\A'hn 



D A/2 



-f 



2 fe A / 2 <S>% / 



(3.9) 



h ^ ) is normalized, so is ^ (see the Appendix [E| for more details). 
We will use the notation: 



Notice that if ^ = 





(3.10) 



where the eigenfunctions ^( Q ), with a £ {h,f}, satisfy the following orthogonality rela- 
tion: 

roc 

S a p5(fi - »') . (3.11) 



^e A M^ (Q) (r)^, (/3) (r) 



Let us therefore turn to the spectrum of the partner operator H.2- Already from the 
pictures of the Schrodinger potentials determining the wave functions (see Fig. [l]) it tran- 
spires that the spectrum of TL2 is, as expected, purely continuous. As for the eigenfunctions, 
even if the partner system ( |3.8] ) is considerably reduced in complexity compared to (^) , 
they cannot be calculated exactly for non-zero \i. It is however possible to find the solu- 
tions of (3.8) in the thin wall limit ( |2.14 ). The determination of solutions being somewhat 
lengthy, we give here only its main points and the results relevant for the forthcoming 
considerations, refering the interested readers to the Appendix [C] for the details of the 
calculations. 

The solution becomes particularly elegant and straightforward for the upper of the 
equations ( |3.8[) , determining the evolution of h. Indeed, in the thin wall limit this equation 
becomes: 

-h" + 8<e 2 siga(r)hf + (I2£ 2 <5(r) - 12£ 4 ) h = fi 2 e 4 ?^h (3.12) 
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Figure 1: Behavior of the Schrodinger potentials of the partner system: to the left potential for h, 
to the right potential for /. The plots are made in the conformal coordinates, w — J r dr ' e~ A{ - r )/ 2 . 
The parameter values are (3 = 1CP 2 , a = 10. 



and its general solution for arbitrary ^ can be found. It can always be written as a linear 
combination of an even and an odd solution. It is the even, appropriately normalized, 
solution which will be of interest for what follows. It is of the following form: 



(3.13) 



with the normalization constant given by: 



2£ 



-1/2 



(3.14) 



m ) * K^e 

The second equation, which determines the evolution of f^, becomes for |r| 3> I /a: 

-/" + 8£ 2 sign(r)/' + (4a 2 + 8£ 2 a - 12£ 4 ) / = ^e^f . (3.15) 

Again, we can easily find the solutions of this equation, however they cannot be extended 
towards r = as simply as the solutions. The (approximate) global solutions of this 
equation can be found for small fx through matching of the \r\ » 1/a solutions obtained in 
the thin wall limit to the exact fi = solutions. In what follows, we will need the odd / M 
(again, properly normalized): 



f [ °\r) 



2£ 



v I ejx 



2vr \Av£ 2 



-3A/2 



\r\ < d 



sign(r)e 4 « 2 Hj^-^ e ^l^ | r | » I , 



(3.16) 



where d ~ 0(£ 2 ). The latter solution is valid for /i < v^ 2 ps a. 
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4. Static potential between sources on the brane 



In the presence of other matter fields, the system of the gravitational and scalar perturba- 
tions couples to the sources via the interaction term: 

1 



S u 



d xdr^fg 



':h MN T MN + 



(4.1) 



Supposing that the coupling of these matter fields to the gravity is much weaker than to 
the scalar field and its perturbations, we restrict our attention to the sources of the field eft. 
This means that our scalar system ( |3.2[) is coupled to the sources of the form J = ( j ) : 



2A'6' + 



2 AY - (2 A" + A ,2 )x + J— (2$" + $'A') = 0, 
— — h 2A — A" 



+ ^-(2^" + ^A') X = J f (x,r) . 

(4.2) 

In this section, we will study the behavior of the potential between two static sources Jm 
and J( 2 ) of the field (j>. 

We suppose that the sources are pointlike and situated at r = r' = 0: 

J(i)(x,r) = J(i) 5 3 (x - xi)(5(r) , 

l(2){x,r) = J(2)5 3 {x - x 2 )5(r) . 

The static potential between Jn\ and J/ 2 ) receives Yukawa-type contributions from all 
modes and is given by the following integral (see Appendix [D] for the calculations leading 
to the expression below): 



V(\x 1 -x 2 \) = - 



dfjL- 







-fJ,\Xl-X2\ 

X\ - x 2 \ 

oo 



E^1)^)( )^W(°U(2) 

a 

-fi\x 1 -x 2 \ 







a) 



, _ ■ , ■ ( 4 - 3 ) 

where ^(a), with a S {^,/} are the eigenfunctions of the operator TCi satisfying the 
normalization condition Q3.11[ ) and the following boundary conditions at r = 0: 

[* Ma) ]l(0) = ^ (a) (0)=0, 
9r[*M«)]3(°)=^(«)( )=°- ( 44 ) 

These are simply 



o 



and 



J m 



(o) 



with /i^ given by ( |3,13j ) and //i 0,1 given by ( 3.16| ). For these we have (see Appendix : 

-1/2 



/,(/)(«) -( 



3 1\ 1/2 /2^ 2 \ 3/2 



7T 2 /3 

Vl6eJ v 4£ 2 f 



f 2 

.7 2 



2£ 2 



+ K 2 



2£ 2 



y-l/2 



1/2 

16e/ V4a 



9z; \ //xe\V/3 



(4.5) 
(4.6) 
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The first expression is valid for an arbitrary p, the second one for p < u^ 2 ~ a. In all this 
interval f^rh is very strongly suppressed and its contribution to the potential can be safely 
neglected compared to the contribution of f^foy Considering distances R 3> 1/a, we can 
keep infinity as the upper limit of integration. The potential therefore becomes: 



T T 3a 



dp 



e-v R 3 1 /2f 



R 7T 2 V H 

p -(2eR)a 



2£< 



2? 



da 



a 3 [J 2 {a) + Y 2 {a) 



(4.7) 



where R = \x\ — x%\ and where we have introduced a dimensionless var iable a = p/(2£ 2 ). 
Remarkably, apart from the factor in front of the integral (identical to the non-gravitating 
case) the behavior of the static potential ( [4,7| ) does not depend on the details of the brane 
structure, but solely on the geometry of the space-time. Indeed, the integral in (fO|) 
contains only one parameter, 2^ 2 R. This means that it is the AdS radius of the five- 
dimensional space-time, (2£)~ 2 , which will determine how the potential V(R) changes 
with the distance between sources. 

Before actually estimating V(R), it is worthwhile to take a closer look at the spectral 
density function p(a) = a~ 3 [./f ( a ) + ^2 2 ( Q )] > which should reveal whether there are any 
resonances present in our system. The behavior of this function is presented on Fig. ||. The 
bulk of the weight of p(a) is concentrated in the region of fairly light modes, p ~ 0(£ 2 ) 
(or a ~ 0(1)), the width of the peak being of the same order as its distance from the 
origin. Thus, from the point of view of the four-dimensional physics the collection of 
modes with masses p ~ 0(£ 2 ) behaves as a wide resonance. This phenomenon is unlike 
any others previously observed in the context of brane scenarios. Indeed, let us remind that 
in the case of quasilocalized massive scalars considered in [23] the spectral density function 
exhibits characteristic features of a narrow resonance: a sharp peak around p = m/y/2, m 
being the mass of the bulk scalar; the width of the peak is suppressed by a factor (to/£ 2 ) 2 . 
The signature of a quasilocalized massless particle would be a density function steeply 
concentrated in the vicinity of the origin [26,27]. 

0.35: 
0.3: 

0.25: 

0.2: 

p(a) 

0.15 : 
0.1: 
0.05: 

5 10 15 20 25 

a 

Figure 2: Behavior of the function p(a) = a~ 3 [Jf ( a ) + ^2 2 ( a )] 
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Let us now proceed to determine how the potential (|4.7|) depends on the distance 
between sources (to simplify the notation, from now on we will set Jmj = J(2).f = !)• 
The dominating behavior of V(R) can be determined at the opposite ends of distance scale, 
for 2£ 2 i? > 1 and for £ 2 /a < 2^ 2 R < 1. At large distances, 2£?R > 1, only small a's 
will contribute, and we can approximate the integrand using the expansion of the Bessel 
functions for small argument. The potential V(R) becomes then: 



V{R) 



3a 
16^R 

3a 
16^R 

3a 
16^R 



a 



2 + a 2 



ci(2V2~fR) cos(2v / 2^ 2 J R) - si(2V2~fR) sin(2\/2£ 2 i?) 



1 1 



2 (2£ 2 R) 2 2 (2£ 2 i?) 4 



+ 



(4.8) 



where the last line was obtained by asymptotic expansion, using 2^ 2 R 3> 1. Hence, in 
this region of distances the static potential between two sources has a 1 /R? fall-off and its 
amplitude is strongly suppressed. The approximate expression for ( |Q| ) is plotted on Fig. || 
together with the numerically evaluated V(R), demonstrating that the analytical result is 
an excellent estimation of the potential's behavior at large distances. 

Let us now estimate the potential at small distances, £ 2 /a <C 2£ 2 i? <C 1. We can 
rewrite (fOI) in the following way: 



V(R) 




da ■ 



a J [Jf (a) + Y%(ot) 



+ 



da 



-2£ 2 Ra 



[,J 2 (a)+Y 2 (a)] 



Q" 



da 



7T 



a*[Jl(a)+Y 2 (a)] 2 



+ 



da 



e -2i 2 Ra _ 1 



TT 



da [ 1 — < 
1 



-2£ 2 Ra 



1 + a 2 



a 2 [,J 2 (a)+Y 2 (a)] 2 a 2 + 1 



2£ 2 Ra 



3 a 
2vr 3 i? 



da 



2£ 2 ,R 



7T 



a3[J 2 2 (a)+y 2 2 (a)] 2 



da 



da ( 1 — e 



a 



1 + a 2 



7T 



[J 2 (a)+y 2 2 (a) 



a< 



2 a 2 + 1 



(4.9) 



While the first integral in (|4.9| ) does not seem to be tabulated, its numerical evaluation 
indicates very suggestively that: 



da ■ 



1 



7T 



a3[J 2 2 («)+y 2 2 (a)] 8 • 
The second integral can be expressed through the sine and the cosine integrals: 



(4.10) 



da L-*e^ _ i 



i 



l + a 2 



7T 

2 L 



7T 



ci(2^ J R) sin(2£ 2 i?) - si(2£ 2 i?) oos^R) - - 



ir$, 2 R\n(£ 2 R) + tt£ 2 R (ln(2) + 7 - 1) 
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Figure 3: Behavior of the potential as a function of R: to the left for the large distances, to the 
right for the short distances. The continuous line is the numerical evaluation of the potential, the 
dotted line represents the analytical estimation. 



where the second line was obtained using the series expansion. Inserting these results 
into (WM) we obtain 



V(R) :! " 



16vri? 



1 + -f RHF R ) + -£ 2R ( 7 " 1 + ln(2) - -Ci 

7T 7T \ 7T 



(4.11) 



where the constant C\ stands for the result of the third integral, which must be found 
numerically: C\ ~ —0.804. The comparison between this expression and the numerically 
evaluated V(R) is presented on Fig. ||. The way in which the above expression depends 
on the strength of gravity is very natural: as £ 2 = /3a, for a fixed a, the smaller is the 
gravitational constant, the larger the interval £ 2 /a <C 2£ 2 R <C 1 and the closer the above 
expression for the potential V{R) becomes to 3a/(167ri?), which is the form of the static 
potential in the zero-gravity (/? = 0) case. When the gravity is switched off, the structure 
of the spectrum changes radically: the fields \ an d (f> get completely decoupled and <f> has 
a normalizable zero mode: 

k{r) = — 2 uq(x) . 

£ cosh (ar) 

The potential mediated by the field uq between two static pointlike sources on the brane 
is in this case: 



AnR lJ y J1 16vri? 

The potential ( |4.7| ) has therefore a continuous limit towards the non-gravitating case 
(see Fig. To conclude, let us observe that the fact that we find the 1/R term in 
the region of distances small compared to the AdS radius is not at all unexpected, as at 
these scales the space-time looks practically flat and the background solution ( 2. 11 )— ( 2. 13| ) 
is very close to the setup of Ref. [1]. 
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Figure 4: Behavior of the potential at short distances: the dotted line is the (3 = case, the 



continuous line represents the gravitating case. The parameter values are (3 — 10 



10. 



5. Periodic source on the brane 

In the previous section we have considered static sources. However, in the actual fact <^'s 
are four-dimensional particles and the source term in ( [4.1] ) is simply an effective description 
of their interactions with other fields localized on the brane. In order to get some insight 
into the dynamics of these particles, we must extend our study to time-dependent sources. 

In this section, we will study the propagation of signal emitted by an oscillating source 
of cf). Let us therefore consider a pointlike source of the form (j f ) situated on the brane 
(more specifically, we consider the hypersurface r = 0) and oscillating with frequency uj. 
The field induced on the brane is given by the convolution of the source with the Green's 
function: 

/oo 
d(t' - t) G(x, r = 0; x',r' = 0) e - *"^ , (5.1) 
-oo 

where G(x,r;x' ,r') is the retarded Green's function of Eqn. (|4.2| ), given by 



d 4 k Mx ^) r,„ V- *M«)( r )*Ma)( r ') 

/i 2 — k 2 — 2ieko 



(5.2) 



where, as in the previous section are the eigenfunctions of the operator Hi, satisfying 

the boundary conditions (^4). As h^(r) components are odd, purely /-type field will 



be produced; it is determined by the (/, /)-component of the Green's function. After 
inserting Grjj\(x, r; x', r') and simplifying we obtain 



r (/J) 



1 r°° pi^^R „ 
(R) = Gft,) (x,r = 0;x',r> = 0) = — d M — - J] [f^ (0)] , (5.3) 

Jo a 



where R = \x — x'\ and where uj^ = \J I J 2 — /i 2 when fi < uj and uj^ = i\J /x 2 — uj 2 
when fj, > uj. Only modes with \x < to are actually radiated; the other ones exponentially 
fall off from the source. Being interested in the propagation of waves, we can integrate the 
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expression for G% f)(R) only up to \x = to. Limiting our interest to frequencies w « a, we 
can neglect the contribution of f^jy given by (f4.6[ ), compared to the one of f^fy, given 
by (|Q|). GVfjJR) becomes then: 

t 2N 3 



G 



, i f)( R ) ~ ^ 



e i^R 3 i 



2 JL 

R vr 2 (3 V A* 



f 2 



2£ 2 



3a 

27T 3 i? 

3a 
2n 3 R 



o 



"/ 2 « 2 da ex P 

~~ 3 



a- 



J 2 (a) + y 2 2 (a) 



da p(a) exp iioR\J\ — (2£ 2 /uj) 2 a 2 



The above expression for G^ ^ (i?) indicates that the behavior of the wave's amplitude will 
depend on its frequency. Let us start by determining this behavior for frequencies uj 3> £ 2 
(or, in other words, wavelengths much smaller than the AdS radius, (2£)~ 2 ). Given that 
the bulk of the weight of p(a) is concentrated in the region of a < O(oo), we can safely 
extend the integration to the infinity and expand the square root in the exponential to 
obtain: 

3a e iulR f°° da ex P 



'(/,/) 



(R) 



-i^a 2 



2vr 3 R J a 3 J 2 (a)+Y 2 (a) ' 



(5.4) 



For distances R <C u/£ , the phase factor varies very slowly and can be set to 1. We obtain 
then: 



G 



(/,/) 



(R) 



3a e 



2?r 3 R 



da 



1 



3a e' 



jR 



a: 



J 2 2 (a) +Y 2 2 (a) IQtt R 



(5.5) 



where we have used ( [4.101) . We therefore recover the usual 1/R dependence of the wave 
amplitude on the distance to the source. 

At very large distances, for R ^> oj/£, 4 , the phase factor varies very rapidly and only 
the small a's will contribute significantly to the integral. We can therefore approximate 
the Bessel functions by the first terms of their series expansion: 



'(/,/) 



(R) 



3a e 



iujR 



167T R 



da 



a 2 + 2 



exp 



.2£ 4 i? 



-a 



to 



By successive integration by parts, we can express this integral as a series in 



^R- 



'UJ) 



(R) 



3a e 



iujR 



16vr R 



UJ 



8^R 



-471-/2 



+ o 



UJ 

~£}R 



(5.6) 



that is, the amplitude is proportional to 1/R 2 , signaling the dissipation of the waves into 
the extra dimension. 

Let us now consider low frequencies: oj <C £ 2 . For these we can simplify the integrand 
using the small argument approximation of the Bessel functions: 



G 



UJ) 



(R) 



3 a 
32^R 

3auj 
128vr£ 4 



-171-/2 



da a exp iojR\/l — (2^ 2 /uj) 2 a 2 

iujR / ,.,T}\ ^iuR/2 



R 2 



2 sin 



ujR \ i 
~2 



ujR 3 
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It is transparent that in this regime of frequencies, the leaking into the fifth dimension is 
important at all distances. Dominating behavior at large distances is again of l/R 2 type. 



6. Conclusions 

We have studied the interactions of the scalar perturbations with sources on the brane, 
both static and oscillating. Our analysis reveals that at distances smaller than the AdS 
radius, R <C £~ 2 , the light modes of the continuous spectrum induce a static potential which 
tends smoothly to the zero gravity limit. Studying the propagation of a signal emitted by 
an oscillating source on the brane, we have shown that for wave frequencies higher than 
the inverse of the AdS radius, the wave amplitude has the usual four-dimensional behavior 
in a large interval of distances, scaling as the inverse of the distance to the source. The 
dissipation of the wave into the fifth dimension becomes significant only at very large 
distances. 

We relate these phenomena to a wide resonance, which can be seen as a remnant of the 
scalar zero mode present in the domain wall setup in the absence of gravity [1] . Presence of 
this resonance ensures a continuous change of the physical quantities when the Planck mass 
is sent to infinity. Provided that the thickness of the wall is much smaller than the AdS 
radius of the space-time, the parameters of this resonance do not depend on details of the 
domain wall's structure, but only on the geometry of the space-time. 

We did not derive here an explicit low energy effective action for our setup. In the 
non gravitating limit, we expect it to coincide with the Nambu-Goto action [28]. The 
translational zero mode can be then identified with the branon excitation. Our results 
show that the gravitational backreaction of the brane on the bulk geometry has a strong 
influence on the appearance of the mode. The effective action for our self-gravitating model 
together with phenomenological implications of our results and their relation with massive 
branons are under study and will be discussed elsewhere. 

Acknowledgments 

We wish to thank S. Dubovsky and A. Neronov for helpful discussions. This work is 
supported in part by Swiss National Science Foundation. 

A. /x = solutions of the eigenvalue equation of r K 1 

In this appendix we present the equations which allow to determine the solutions of the 
eigenvalue equation of the operator TCi for fj, = and its general solution. In this case, the 
general solution of ([3^) can in fact be determined exactly. 

First of all, as suggested by the factorization of Hi, we get two solutions solving: 



A 




,3A/4 [ n 




) 



= . 
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The above system of first order equations can be easily decoupled by eliminating one of 
the variables, giving e.g.: 

h' ' + [A 1 - 2— J h' + \2A" - 2A'— J h = 0, 



whose solutions are: 









l-( 




l-( 


U 2 V 





A'e 




— $'e A 
3 



1- A'e 



dye A ^ 



The remaining two solutions can be found by solving an inhomogenous equation 

A 



h 
Jo. 



h 
Jo y 



(A.l) 



where the inhomogenous term is the solution of 



0, given by: 
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Jo 
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^ -3A/2A' 

2k $" 



-3A/2 



dye J 
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-3A/2 



.4' 



The two additional solutions coming from ( A.l ) are of the following form: 
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dr\e 
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+ 



+A'e~ A 



\ V 2k 
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n dr 2 e- 2A ^ [2$" + AV] f'drse 



+ 2 / dne J 



L dr 2 e- 2A ^ [2$" + A'$'] 



None of these solutions is normalizable. 
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B. Connection between the spectra of Tix and 7t 2 



As long as f3 7^ 0, the operators TL\ and TL2 have both purely continuous spectra. In this 
appendix we will illustrate in some detail how the eigenfunctions of the two operators are 
related. 

Suppose that we have an orthogonal basis of the eigenfunctions of the operator Ti.2 
for some boundary conditions. As the partner system is decoupled, the natural basis is 
composed of 



n(h) 



hp 




and \& 



where, trivially, 'j? ^-type solution vectors are orthogonal to ^ / A-type ones. The basis 
of the eigenfunctions for the coupled system can be obtained from these using ( |3.9[) : 



n(h) 



1 

H 



h, 



and * 



' 



where again )]/ ,^%-type solution vectors are orthogonal to Vl/^^-type ones. Indeed, we have 
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00 r n T r 

are ^ W M(h) y\. w M#(/) 

-00 



/i- 



4t / rfre A W AA + * 



MM 



dre A(r)^T ^ =0 . 



Moreover, the following relations are satisfied: 



Kf) /"(/) 



K/) /"'(/) 



/oo 

00 

Consequently, normalizing the basis solutions of the partner system using: 

dre A (%(r)V(r) =*(//-//') , 

dre A «/ M (r)/^r) =%-//) , 
we obtain an orthonormal basis of eigenfunctions Wi: 

^e A ( r %J (Q) ^, (/3) =6{fi-fM')S aP , 

where a, (3 G {/t, /}. The solutions ^Va) satisfy the following completeness relation: 



00 
00 



/>oo 

5^ / dM*„ (a) (r)tf 
^ ^ 



2 • 



(B.l) 
(B.2) 

(B.3) 
(B.4) 
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C. Eigenfunctions of the partner operator 7i 2 

This appendix contains a detailed derivation of the eigenfunctions of the partner operator 
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For /j, = the general solution of the system Q3.8Q can be calculated exactly: 
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There are two combinations of the four basis solutions which are solutions of A* ( \° ) = 



,(0) 



fo 



D 



(o) 



These correspond to D 3 — y q— emu. ^ 4 — — y 2^-^ 2 

For non-zero [i it is impossible to find exact solutions. However, one can find approxi- 
mate solutions of the system using the thin wall limit. The partner system being decoupled, 
our task is reduced to finding solutions of two second order equations. 
• Let us start by determining the solutions of the first of the equations (|3.8| ) : 



-h" - 2A'h' 



A" + ^A' 2 )h = f i 2 e- A h. 



(C.l) 



The exponent of the warp factor A(r) is given by Eqn. fl2,12| ). In the thin wall limit (2.14), 
we have 



A 
A' 
A" 



< 2 \r\ , 
-4£ 2 sign(r) 



Consequently, Eqn. (|C.l[) becomes : 

-h" + 8f siga(r)h' + (I2e 2 5(r) - 12£ 4 ) h = ^e^h , 
where £ 2 = a(3. This equation can be solved for arbitrary fi and its solutions are: 
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Eqn. being symmetric under reflection r — > — r, for any given \x there exists an even 

and an odd solution. Since they are, of course, linearly independent, the general solution 
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can be written as a linear combination of them. We can therefore consider solutions with 
definite parity. For an even solution, we must have B3 = B\ and -B4 = B2 and it can thus 
be written as: 



r = e 



4e\r\ 



2C 2 



(C.3) 



For h^ e \r) to be a solution of ( |C.2| ), it is necessary that the coefficients Bj 6 and B, 
satisfy the following relation (imposed by the jump of the derivative at r = 0): 
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(C.4) 



Imposing the normalization condition ( |B.1[) , we fix the remaining constant: 



4£2 



1 + 



(C.5) 



For an odd solution, we must have U3 = —B\ and 1?4 = — B^. Thus, it can be written as: 



U°Xr) = sign(r)e 452|r| 



(C.6) 



The coefficients and .B^ must satisfy the following relation (imposed to avoid <5'-like 
singularities) : 
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Imposing the normalization condition flB.l|), we find: 
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(C.8) 



Let us now consider the second of the equations (^ 
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In the thin wall limit, for |r| 3> 1/a, this equation becomes: 

-/" + 8£ 2 sign(r)/' + (4a 2 + 8£ 2 a - 12£ 4 ) / = /.V^l/ , 
the solutions of which are: 
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D 3 J U 



2? 



+ D 4 J. 



2e 



r > , 
r < , 



(C.ll) 
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where v = l/f3 + 1. We cannot extend them towards r = as simply as it could be done 
for the solutions of Eqn. ( |C.2| ). However, it is still possible to find (approximate) global 
solutions of ( C.10| ) for small /i matching ( p,ll| ) to the exact \i = solutions. 

Again, given the reflection symmetry of ( p.9| ) we can consider even and odd solutions 
separately. We will limit ourselves to determining the odd solutions, as only these are 
needed in the calculations performed in sections |3| and || . The even solutions can be found 
in an entirely analogous manner. The odd /j, = solution is: 



JO 



(o) 



-3A/2 



A' 



For small /Lt, the global odd solution of ( |C,10 ) can therefore be approximated by: 



C (o) e -3A/2^_ 



sign(r) e 



4£ M 



\r\ < d , 
1 

r > - 
a 



where d must be such that the /x 2 (and higher) corrections to the /i = term are negligible. 
We must have d < (2f)~ 1 ln(2£ 2 ufj, ). In order to be able to do the matching, we can 
safely take d ~ 0(£ -2 ). Then our solution for is valid for fi < v^ 2 ~ a. For |r| S> 1/a 
the exact solution can be expanded in series in e _2a ' r ': 



/(°) = e" 3A / 2 ^ « - sign(r)^ e 2 «M+ 6 « 2 M (l + 2 e - 2 ^ + 0{e~^ 



f 2 

-sign(r)- 

and in the overlap region a -1 <C \r\ < d we have 



,4€ 2 |r| x / ^2e\r\\ „ _1 



" ' 2£ 2 



e ii 

2ttv V^i 2 

1 / e ;U 



2-kv \Au^ 



a 2a|r|+6£ 2 |r| 



,-2a|r|+2£ 2 |r| 



Comparing the above three expressions, we obtain (at the lowest order in [i) D^=0 and 



c (o) 



va 



e fi 



2ttv£ 2 \4v£ 



2tt \4^ z 



(o) 



Therefore, in the first approximation, the solution is 



/»(r)«Z?r< 



to 



efi 

2vr V^K 2 



-3A/2^ 



Irl < d , 



sign(r)e 4 « 2 Hj^-^ e 2 « 2 M^ | r | » I 



(C.12) 



Imposing the normalization condition ( p.2[ ) we fix the remaining constant 

2 



to. 



4^2 



(C.13) 



which completes the determination of f(°\r 
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D. Determination of the static potential 

In this appendix, we present the details of the calculation of the potential between two 
static sources of (j). The field produced by the source J is given by: 

^r)=/^Q(„r,^, (D.l) 

where G(x,r;x',r') is the (matrix) retarded Green's function of the coupled system (|4,2|), 
satisfying 



e-Vd/Aj 



G(x, r; x', r') = 5 (4) (x - x')<5(r - r') e - 2A( - r) 1 2 . 



The Green's function G(x,r;x',r') can be written as a Fourier transform of the Green's 
function G^(r, r') of our one-dimensional matrix Sturm-Liouville problem ( |3.2[ ): 

(Wi-A ; 2 l 2 )Gf 5) (r,r / )=^(r-r / )e- A Wl 2 . (D.2) 

As per usual, G/U (r, r') can be expanded in the basis of eigenfunctions a E {/i, /} of 

the operator Hi, which allows us to write the following expression for the retarded Green's 
function: 

G(x,r;x',r') = j ^ e ^-^) G f 5) (r, r') 

J (2vr) 4 Jo * ^ ^_ k 2_ 2i£ k y 1 

If the source J is static, we can perform in ( p.l ) the integration over time on the Green's 
function and obtain 

#(x,r) = y d 3 x'dr'e 2A ^ ^J°° dt'G{x, r; x', r') \ J(f, r') 

1 /* /i ^ | 2* JJJ ' I 



The static potential is the interaction energy of a source = ( jW J with the field ^(2) 





produced by J( 2 ) = j( 2 ) ) , that is: 



y = y d 3 xdre 2A( - r) f[ 1) {x,r)^ i2) (x,r) 

d 3 xdre 2A ^ J d 4 x'dr'e 2A ^jJ 1) (x,r)G{x,r-,x',r')J {2) (x',r') 



— I d 3 xdre 2A ^ I d 3 x'dr'e 2A( ^x 
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roo e -fi\x—x'\ 

j ^ E4)( i ' r )^)( |, ) $ i«)(^)(^/ 
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Suppose that the sources are pointlike and situated at r = r' = 0, that is 

I{i)(x,r) = J ( i) <5 3 (x - £i)(5(r) , 
J(2)(x,r) = !{2)^{x - x 2 )S(r) . 

Then, the static potential V becomes: 



V(\£!-X2\) 



1 



-fi\Xl— X2\ 



4vr J Q \xi-x 2 \ 



2) 



1 „ 



e -fi\x 1 -x 2 \ 



where in the last line we have used the fact that both sources are the type ( j f 
E. Solutions of the coupled system 

In this appendix we determine the solutions of the coupled system Hi, using its connection 
with its (decoupled) partner system H 2 via the relation ( |3.9| ). 

Applying A* on the even solution of the eqn. ( C.l ), given by ( p,3| )-( |CT5| ), we obtain: 
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We have therefore: 
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where we have used (|C.5| ). 

Applying A^ on the odd solution of the eqn. ( C.9 ), given by ( C.12|) -( [C.13[ ) we obtain: 
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which for \r\ < d gives us: 



6P ( en \ v D<f\v) 



A'e 



i„—A 



/i 



6/3 /We 



while for r > 1/a we have: 



The expression ( |E.1| ) gives us: 



(E.l) 



v I efj, 



16e V4^ 2 ^ 



2v-l 



where we have used (C.13). 
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